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In this work, we numerically demonstrate an all-optical tunable Fano resonance in a fishnet 
metamaterial(MM) based on a metal/phase-change material(PCM)/metal multilayer. We show that the 
displacement of the elliptical nanoholes from their centers can split the single Fano resonance (FR) into a 
double FR, exhibiting higher quality factors. The tri-layer fishnet MMs with broken symmetry accomplishes 
a wide tuning range in the mid-infrared(M-IR) regime by switching between the amorphous and crystalline 
states of the PCM (Ge 2 Sb 2 Te 5 ). A photothermal model is used to study the temporal variation of the 
temperature of the Ge 2 Sb 2 Te 5 film to show the potential for switching the phase of Ge 2 Sb 2 Te 5 by optical 
heating. Generation of the tunable double FR in this asymmetric structure presents clear advantages as it 
possesses a fast tuning time of 0.36 ns, a low pump light intensity of 9.6 fiW/fim 2 , and a large tunable 
wavelength range between 2124 nm and 3028 nm. The optically fast tuning of double FRs using phase 
change metamaterials(PCMMs) may have potential applications in active multiple-wavelength nanodevices 
in the M-IR region. 



Recently, Fano resonance(FR) excited by an interference coupling between bright(radiating) modes and dark 
(non radiating) modes in metallic plasmonic nanostructures has been widely investigated 1,2 . It has gener- 
ated much interest because of its application in the fields of chemical and biological sensing 3 ' 4 , surface 
enhanced raman spectroscopy 5 , active plasmonic switching 6 , and slow light devices 7 . FR exhibits an asymmetric 
line shape with a sharp feature in transmission, reflection or absorption spectrum over an extremely narrow 
frequency range 8 ' 9 . There have been a few attempts to create FRs in metallic plasmonic nanostructures, such as 
core-shell structure 10 , metallic nanowire arrays 1112 , composite cut-wire structures 13 plasmonic dimmers/clusters 14 
and photonic metamaterials(MMs) 15 " 18 . Among the various MMs, the configurations with broken symmetry are 
the most common approach for exciting narrow FRs. These asymmetric MMs give rise to a non-uniform 
electromagnetic environment around the resonators, leading to the effective coupling between bright modes 
and dark modes 19 . 

The ability to control the central wavelength of a FR in MMs can find more flexible and critical applications in 
practice such as tunability of dispersion and group velocity of light 18 , multispectral imaging/detection 20 ,broad- 
band delay lines 21 . However, most researchers focus on passive tuning of FRs by adjusting geometric parameters of 
the MMs 22 " 26 which limits the scope of applications. Recently, active tuning of the optical response of MMs has 
attracted growing interest owing to its ultra-short modulation time and on-to-off switching capabilities 27 ' 28 . 
Therefore, for the purpose of actively controlling FRs, MMs integrated with optically active materials is a 
promising approach. Quite recently, a variety of approaches to modulate FRs in active MMs have emerged. 
Pryce et al exploited tunable FRs in MMs by means of mechanical deformation of an elastomeric substrate 29 . 
Belotelov et al presented a magnetically tunable FR in a magneto-plasmonic MM 30 . Kurter et al reported cryogenic 
temperature-dependent electromagnetically induced transparency(EIT) MM 31 Samson et al demonstrated a 
frequency shift of a FR mode in a plasmonic planar MM by electrically tuning the dielectric properties of an 
adjacent chalcogenide glass layer 32 . Mousavi et al used the plasmonic response of graphene to tune a FR in the 
mid-infrared(M-IR) region 33 . However, FRs in these MMs are magnetically and electrically tuned hence only 
responding on a scale of several microseconds to seconds 34 . In addition, the integration of the required MM 
structures as well as the electrodes for the tuning of active dielectric layers can be difficult. Therefore, an appealing 
route that is notable for fast tunable FRs is to optically adjust the optical properties of the MMs. To this end, Gu 
et al presented active optical control of EIT MMs by integrating photoactive silicon(Si) islands into resonant 
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elements 18 . However, this MM works only in the terahertz(THz) 
regime, and its structural geometry is complicated to fabricate in 
the optical region, which limits its potential use. Zhu et al proposed 
a fast all-optical tunable FR in the optical frequency range, using a 
plasmonic MM made of gold and polycrystalline lithium niobate 
(LiNb0 3 ) 8 . However, this two dimensional (2D) planar MM system 
is based on a single metallic grating layer. It has been suggested that 
multiple layers of perforated metal-dielectric-metal(MDM) stacks 
with up to 100 to 200 layers along the propagation direction show 
a promising approach for three dimensional (3D) optical MMs 35-37 . 
Such 3D MMs can possess a thickness much larger than the free 
space wavelength in the optical region and excite a strong mag- 
neto-inductive coupling between neighboring functional layers 
under a normally incident light. The tight coupling between adjacent 
LC resonators through mutual inductance results in a low loss 36 . 
Thus 3D multilayer MMs are likely to be required for nearly all 
potential applications, such as cloaking and invisible materials 35-38 
as well as tunable metamaterials. Here we choose phase change 
material(PCM) as the dielectric interlayer,using which,we have prev- 
iously shown designs for low loss tunable negative index material 39 
and a polarization independent tunable perfect absorber 40 .Thus, 
PCM based active tuning of FRs in three dimensional(3D) multilayer 
MMs is an interesting direction to pursue. To the best of our know- 
ledge, it is still a formidable challenge to obtain a fast, low-power all 
optical tunable FR using 3D multilayer MMs in the optical regime. 

Here, we numerically demonstrate that the resonant frequency of 
FRs can be optically tuned using PCMs in the mid-infrared(M-IR) 
region. Our structure is composed of an elliptical nanohole array 
(ENA) perforating through a Metal-Dielectric-Metal(MDM) tri- 
layer, and a prototypical PCM: Ge 2 Sb 2 Te 5 , is selected as the dielectric 
layer. A single FR is observed in the symmetric MDM-ENA. This FR 
response is caused by the interaction between the bright and dark 
modes. The elliptical nanohole resonator array shows electric res- 
onance (namely localized surface plasmon) which excite the bright 



modes, whereas MDM multilayers sustain a magnetic resonance 
(namely inductive- capacitive resonance) to induce the dark modes 
at the same frequency 18,41 . As soon as the elliptical holes are displaced 
from their central positions, the single magnetic resonance would be 
divided into two magnetic resonances at different wavelengths. A 
double FR appears, owing to the destructive interference between 
the electric resonance and the two magnetic resonances at different 
wavelengths. This interference can enhance the FR. By using 
Ge 2 Sb 2 Te 5 , a broad spectral tunability (from 2124 nm to 3082 nm) 
of the double FRs in the M-IR region can be obtained by switching 
between the amorphous and crystalline states of the PCMs 42 . 
Importantly, a heat model is constructed to examine the temporal 
variation of the temperature of the Ge 2 Sb 2 Te 5 layer in the asymmet- 
ric MDM-ENA. The model shows that the temperature of the 
amorphous Ge 2 Sb 2 Te 5 layer can be raised from room temperature 
to >433 K(phase transition point of Ge 2 Sb 2 Te 5 ) 43,44 in just 0.36 ns, 
with a low operating pump light intensity of 9.6 fiW/fim 2 . Therefore, 
our proposed structure can rapidly tune the double FR under a low 
light intensity in the optical region using Ge 2 Sb 2 Te 5 . The optically 
active tuning of double FRs in these phase change metama- 
terials(PCMMs) may open up avenues for designing tunable multi- 
wavelength nanostructures in the optical regime such as slow light 
devices, ultrasensitive sensors and nonlinear devices. These PCMMs 
possess a simple geometry which remains compatible with standard 
photolithography patterning and can be easily experimentally rea- 
lized in the optical region. Finally, it should be noted that PCMs do 
not require any energy to maintain the structural state of the mater- 
ial. Thus, once the MM has been switched it will retain its FRs until it 
is switched again. This obviously makes tunable FRs in PCMMs 
interesting from a 'green technology' perspective. 

Results 

Metamaterials design. The proposed MM structures consist of two 
metallic films (30 nm thick Au) separated by a dielectric interlayer 
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Figure 1 | (a) Schematic of the symmetric MDM-ENA((5 = 0 nm) consisting of a 160 nm thick Ge 2 Sb 2 Te 5 dielectric layer between two 30 nm thick Au 
films perforated with a square array of elliptical holes suspended in a vacuum, (b) Illustration of the symmetric element of ENA, the lattice constants are L x 
= 500 nm, L y = 400 nm and hole diameters are dj = 320 nm, d 2 = 150 nm. (c) Schematic of the asymmetric MDM-ENA(£ = 70 nm) consisting of a 
160 nm thick Ge 2 Sb 2 Te 5 dielectric layer between two 30 nm thick Au films perforated with ENA suspended in a vacuum, (d) Illustration of the 



asymmetric element of ENA, the lattice constant is L x 
upper elliptical hole from the centred = 70 nm). 



500 nm, L y = 400 nm, hole diameters are di = 320 nm, d 2 — 150 nm, 3 is the distance of the 
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Figure 2 | Dielectric constant (a) £\(co) vs wavelength, (b) s 2 (co) vs wavelength for both amorphous and crystalline phases of Ge 2 Sb 2 Te 5 . 



(160 nm thick Ge 2 Sb 2 Te 5 ). Here, we simulate two sets of multilayer 
MMs starting from a symmetric element with both top and bottom 
elliptical holes exactly at the centre (see Fig. l(a)-(b)) and then 
displaced the upper elliptical hole from the centre by a distance 
"S" while keeping the lower elliptical hole fixed at the centre in 
order to break the symmetry (see Fig. 1(c) -(d)). For both 
structures, the elliptical nanohole array (ENA) (lattice constant, 
along the x axis is L x = 500 nm, along the y axis is L y = 400 nm; 
hole diameters, d x = 320 nm, d 2 = 150 nm) is perforated through 
the entire Au/Ge 2 Sb 2 Te 5 /Au structure, /? is a cross-section plane of 
the structure, and the elliptical holes are periodically arranged in both 
the x and)/ directions. The Au bottom layer interacts with the upper 
Au layer to give rise to a closed circle of displacement current (J D ) and 
localize an electromagnetic (EM) field within the dielectric interlayer. 
The structures are consider to be suspended in a vacuum, deep 
etching of a silicon support substrate can be used to achieve this. 
Au is selected as the metal due to its stability and low ohmic loss. The 



geometry of the unit cell pattern and the thickness of the sandwich 
layers have been chosen to allow for the impedance matching 
between the metamaterials and impinging plane wave 45 . The 
simulation is performed by commercial software Comsol, which is 
based on the Finite Element Method (FEM). The dielectric properties 
of Au as given by Johnson & Christy are used 46 . All the structures are 
excited by a source with a wavelength range from 1000 to 4000 nm, 
propagating along the negative z direction with the E field polarized 
along the small axis of the elliptical hole as shown in Fig. 1(a) and(c). 
The light source has a recurrence rate, f r = 25 kHz and pulse 
duration of 0.25 ns. The light fluence illuminating the sample from 
a single pulse is written as 47 

F,(r)=-^exp(-^) (1) 

where P 0 = 3 mW is the total incident power, r is the distance from 
the beam center, w = 10 /nm is Gaussian beam waist. To account for 
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Figure 3 | The comparison of (a) the transmittance, (b) the absorbance between asymmetric MDM-ENA((5 = 70 nm) with amorphous Ge 2 Sb 2 Te 5 and 
symmetric MDM-ENA(£ = 0 nm) with amorphous Ge 2 Sb 2 Te 5 . (c) Zoom in picture of Fig. 3(a). (d) Zoom in picture of Fig. 3(b). 
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Figure 4 | The comparison of (a) the transmittance, (b) the absorbance between the asymmetric MDM-ENA(5 = 70 nm) with amorphous Ge 2 Sb 2 Te 5 
and crystalline Ge 2 Sb 2 Te 5 . 



the periodic nature of the metamaterial, the model boundary at 

L x 

x= + — and y = ±L y is set to condition of perfect magnetic 

conductor and perfect electric conductor,respectively for normal 
incidence. Scattering boundary conditions are applied in the z 
direction. Adaptive meshes are generated by a built-in auto 
meshing algorithm of Comsol. The frequency dependent reflection 
r(co) and transmission t(co) are obtained from the S-parameters in the 
simulation package, and the absorbance is calculated as 1 — | r(co)\ 2 — 

Optical properties of Ge 2 Sb 2 Te 5 . The real, 8i(co) and imaginary, 
s 2 ((jo) parts of the dielectric function of Ge 2 Sb 2 Te 5 in different 
phases are obtained from the experimental data in 42 , and for the 
mid-infrared (M-IR) spectral range the dielectric function is 
shown in Fig. 2; as can be seen there is a very large change in the 
dielectric function's real component. The dielectric constant of 
Ge 2 Sb 2 Te 5 is very dispersive and has a non- negligible imaginary 
part. It also changes considerably during the reversible structural 
transformation from amorphous to crystalline. Different PCMs 
can display similar behavior in other parts of the spectrum. It 
should be mentioned that the reversible phase transition in 
Ge 2 Sb 2 Te 5 is highly repeatable and more than a billion cycles have 
been experimentally demonstrated in data storage devices 49 . 

Tunability of double Fano resonance in the phase-change 
metamaterial. Figure 3 presents transmittance and absorbance for 
the structures of symmetric MDM-ENA(c) = 0 nm) and asymmetric 
MDM-ENA(c) = 70 nm) at normal incidence, where the dielectric 
layer is amorphous Ge 2 Sb 2 Te 5 .Here, we show that the symmetric 
MDM -EN A exhibits a single FR (denoted as PI mode) with 
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Figure 5 | 3D-FEM simulation of (a) total electric field intensity distribution 
plane for the amorphous symmetric MDM-ENA((5 = 0 nm) at PI mode(/l = 



resonance wavelength centered at about 2098 nm (898.5 THz)(red 
dash curve in Fig. 3(a) with zoom-in in Fig. 3(c)). Spectral width of 
the single FR is 13.4 THz, leading to a quality factor(Q) of 67, where 
Q factor is defined by the resonant frequency over the width of the 
resonance 50 . It can be inferred that this single FR must be excited by 
the interaction between the electric resonance and the magnetic 
resonance in the structure. By displacing the upper elliptical hole 
from its centre by a distance S = 70 nm, the single FR is split in 
two distinct resonant peaks at 2124 nm(887.5 THz) (marked as P2 
mode) and 2050 nm (marked as P4 mode) by an asymmetric 
resonant dip. Namely, the single magnetic resonance at the 
wavelength of 2098 nm would be divided into two magnetic 
resonances at different resonance wavelengths of 2050 nm and 
2124 nm by breaking the structural symmetry. Therefore, the 
double FR in the transmittance spectrum(black solid curve in 
Fig. 3(a)) is attributed to the destructive interference of the electric 
resonance(bright modes) with the magnetic resonances(dark modes) 
at different resonant wavelengths. In particular, a narrow transpa- 
rency sub -band resonance window with a spectral width of 10.3 THz 
is achieved at P2 mode for the asymmetric MDM-ENA. It gives rise 
to a higher Q factor of 86 than the Q factor of 67 in the symmetric 
structure. This shows that asymmetric elements can significantly 
improve the FR in 3D multilayer MMs. In Figs. 3(b) and 3(d), one 
can observe that the absorbance in the asymmetric MDM-ENA is 
also split accordingly, whereas possessing a similar maximum value 
with the symmetric MDM-ENA. 

Figure 4 shows the evolution of the double FRs in both the trans- 
mittance and absorbance spectrums of the asymmetric MDM- 
ENA^ = 70 nm) for M-IR radiation at normal incidence, in the 
amorphous and metastable cubic crystalline states. With the phase 
change of Ge 2 Sb 2 Te 5 , the double FR at 2124 nm (P2 mode) and 
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Figure 6 | 3D-FEM simulation of (a) total electric field intensity distribution, (b) total magnetic field intensity distribution and J D distribution 
along P plane, at P2 mode(^ = 2124 nm); (c) total electric field intensity distribution, (d) total magnetic field intensity distribution and J D distribution 
along p plane, at P4 mode(A, = 2050 nm) for the amorphous asymmetric MDM-ENA(5 = 70 nm). 



2050 nm(P4 mode) for the amorphous state shifts to a double FR at 
2932 nm (P5 mode) and 3028 nm(P6 mode) for the crystalline 
state. Figure 4(b) presents the trapped plasmon mode absorbance 
lines of the structure for different states of Ge 2 Sb 2 Te 5 .The 
Ge 2 Sb 2 Te 5 variable dielectric function gives rise to a concomitant 



tunability in the metamaterial absorbance. We also find that the 
peaks of both transmittance and absorbance drop down when 
Ge 2 Sb 2 Te 5 changes from the amorphous to crystalline in addition 
to a shift, owing to a reduction in the impedance matching of the 
structure to vacuum. 
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Figure 7 | 3D-FEM simulation of (a) total electric field intensity distribution, (b) total magnetic field intensity distribution and J D distribution 
along p plane, at P5 mode(^ = 2932 nm); (c) total electric field intensity distribution, (d) total magnetic field intensity distribution and J D distribution 
along p plane, at P6 mode(X = 3028 nm) for the crystalline asymmetric MDM-ENA(5 = 70 nm). 
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Figure 8 | (a) 3D -FEM simulation of the transmittance spectrums (red solid) of the symmetric MDM-ENA with amorphous Ge 2 Sb 2 Te 5 and best-fits to 
the FRlineshape Eq.(l) (blue dash); (b) 3D-FEM simulation of the transmittance spectrums (red solid) of the asymmetric MDM- ENA with amorphous 
Ge 2 Sb 2 Te 5 and best-fits to the FR lineshape Eq.(l) (blue dash); (c) 3D-FEM simulation of the transmittance spectrums (red solid) of the asymmetric 
MDM- ENA with crystalline Ge 2 Sb 2 Te 5 and best-fits to the FR lineshape Eq.(l) (blue dash). 



Figure 5 shows total electric field intensity distributions 
E = yj \E X \ 2 + \E y \ 2 + \E Z \ 2 , total magnetic field intensity distribu- 
tions H= y/\H x \ 2 + \H y \ 2 + \H z \ 2 and J D along the ft plane at the 
resonant wavelengths of 2098 nm(Pl mode) for the symmetric 
MDM -EN A( 5 = 0 nm) with amorphous Ge 2 Sb 2 Te 5 . Here, the 
arrows represent J D whereas the color represents the magnitude of 
the electromagnetic field intensities. In Fig. 5(b), the formation of the 
closed J D loops are observed to support a magnetic resonance at 
which light is trapped and strongly absorbed, thus the H field can 
be efficiently confined between two Au layers. However, Fig. 5(a) 
shows a low concentration of the E field intensity within the dielec- 
tric interlayer and elliptical apertures, indicating a weak electric 
resonance. 

Figure 6 demonstrates E, H and J D along the ft plane at the res- 
onant wavelengths of 2124 nm (P2 mode) and 2050 nm(P4 mode) 
for the asymmetric MDM-ENA(c) = 70 nm) with amorphous 
Ge 2 Sb 2 Te 5 . Figure 6(b) shows that H field intensity distributions 
are still confined well within the dielectric layer at P2 mode. 



Table 1 | Parameters used in the FR lineshape to fit the transmit- 
tances of the MDM-ENA 
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Whereas, Figure 6 (d) shows that H field intensity at P4 mode is 
weaker hence leading to a smaller magnetic resonance dipole than 
both PI in Fig. 5(b) and P2 mode in Fig. 6 (b). Figure 6(a) and (c) 
show the E field intensity distributions at P2 and P4 mode, respect- 
ively. As can be seen, the E field intensity distributions at P2 mode in 
Fig. 6(a) are more efficiently concentrated thus giving rise to much 
stronger electric resonances than PI in Fig. 5(a) and P4 mode in 
Fig. 6(c). Therefore, for the asymmetric MDM-ENA, the interaction 
of the stronger electric resonance and magnetic resonance leads to 
the narrower FR at P2 mode compared with PI mode in the sym- 
metric MDM-ENA. 

In Fig. 7(a)-(d),we present the distributions of E, H and J D assoc- 
iated with P5 (X = 2932 nm) and P6 (X = 3028 nm) modes for the 
asymmetric MDM -EN A( 5 = 70 nm) with crystalline Ge 2 Sb 2 Te 5 . It 
demonstrates that the closed J D loops are still maintained at two 
resonance modes. It provides two magnetic dipolar moments thus 
leading to double resonant peaks. However, the double FRs is broa- 
dened owing to the inference between the weak magnetic resonances 



Table 2 | Material thermal properties used in the Heat transfer 
model 

Special heat 

capacity Cs Density thermal conductivity 

(J/(kg-K)) p(kg/m3) k(W/(rrrk)) 

gold 129 56 19300 56 31 /(bulk) 56 , 91 

(thickness = 30 nm) 56 , 
Ge2Sb2Te5 220 55 6150 55 Temperature dependence 55 

air l 47 353 [K] 47 0.03 47 

T 
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Figure 9 | 3D- FEM simulation of heat power irradiating on an 
amorphous asymmetric MDM-ENA located at the beam center, where 
the solid red line presents the heat power irradiating on the structures 
under normal incident intensity of 9.6 /*W//mi 2 ,the dash red line is the 
temperature of the amorphous Ge 2 Sb 2 Te 5 layer during one pulse. 

and electric resonances, resulting from the low intensities distribu- 
tions of H and E fields in the structure. 

Fano formula fitting to the transmittance spectrums simulated 
through FEM method. In Fig. 8, we fit the transmittance curves of 
both the symmetric MDM-ENA(c) = 0 nm) and asymmetric MDM- 
ENA(5 = 70 nm) using the FR lineshape in Eq.(2) 51 " 53 : 

\2 



l + sj 



,('=1,2, 



(2) 



where s z = 



2(co — a>i) 
Ti ' 



particularly, q A are the phenomenological 



shape parameters (so-called asymmetry parameters), coi are the 
resonance frequencies and r b are the widths of the autoionized 
states (the resonance linewidth at half maximum, FWHM) of the 
resonant modes Pi 53,54 . With the parameters shown in Table 1, the 
calculated transmittances based on Eq.(l) approximately reproduce 
our FEM simulation data for the symmetric MDM-ENA(5 = 0 nm) 
with amorphous Ge 2 Sb 2 Te 5 , asymmetric MDM-ENA(5 = 70 nm) 
with amorphous and crystalline Ge 2 Sb 2 Te 5 , respectively. 

It is found that the Fano formula mostly matches the FEM simu- 
lation data around the resonant modes. Q factor is defined by the 
resonant frequency over the width of the autoionized state 15,51 , where 
the resonant frequency and the width of the autoionized state are 
obtained from the Fano formula fitting transmittance curves shown 
in Fig. 8. 



Photothermal effect of the phase-change metamaterial. Since the 
reversible amorphous - crystalline phase transition of Ge 2 Sb 2 Te 5 can 
be induced through optical heating, it is essential to understand the 
heat induced switching behavior of the asymmetric MDM-ENA. To 
show this, a heat transfer model is used to investigate the temporal 
variation of temperature of Ge 2 Sb 2 Te 5 layer at normal incidence 
using the Finite Element Method(FEM) solver within COMSOL. 
The material thermal properties used for the simulation are 
summarized in Table 2. The thermal conductivity of Ge 2 Sb 2 Te 5 
changes with the temperature is obtained from the experiment 
data in 55 . 

In this simulation, the thermal energy absorbed by one unit cell is 
defined as 47 



E th (r)=R a xL x xL y xFi(r) 



(3) 



where L x = 500 nm is the lattice constant of the structure along the x 
axis; L y = 400 nm is the lattice constant of the structure along the y 
axis; R a = 0.0438 is the absorption coefficient of the structure, 
derived from the overlap integral between the light source power 
density spectrum and the absorbance spectrum, shown in Fig. 4(b). 
The heat source power is then described by a Gaussian pulse function 



Q s (r,t) = E th (r) -L- exp ( - <LJ«L 

V71T T z 



) 



(4) 



where t = 0.15 ns is the time constant of the light pulse, t 0 = 0.3 ns is 
the time delay of the pulse peak. Figure 9 shows Q s (r, t) and the 
temperature of the amorphous Ge 2 Sb 2 Te 5 layer, where the structure 
is located at the center of Gaussian beam. The numerical simulation 
shows that the temperature within the amorphous Ge 2 Sb 2 Te 5 dielec- 
tric layer is a function of the incident radiation flux and can exceed 
the amorphous to crystalline phase transition temperature of 433 K 
after 0.36 ns and a maximum of 472 K after 0.59 ns under a thresh- 
old incident flux of 9.6 /aWI/nm 2 . The temperature starts dropping 
after 0.59 ns before the next pulse arrives, due to heat dissipation to 
the surroundings. 

The temperature distributions of the structure at 0.36 ns and 
0.59 ns along the plane /?are shown in Fig. 10(a) and (b) respectively. 
One can observe that the temperature within the amorphous 
Ge 2 Sb 2 Te 5 layer is uniform, and the dominant temperature gradient 
is towards the top and bottom Au films. 

Discussion 

In conclusion, we have numerically demonstrated a low-power and 
rapid all -optical tunable double FR in an asymmetric fishnet MMs 
based on a Au/Ge 2 Sb 2 Te 5 /Au multilayer. By breaking the symmetry 
in the ENA element, the single FR from the symmetric MDM-ENA 
splits into two resonant peaks, resulting in a double FR. This inter- 
ference leads to a sharp transparency peak with higher Q factor. A 
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Figure 10 | The temperature distribution of the unit cell of the amorphous asymmetric MDM-ENA along the P plane at (a) 0.36 ns and (b) 0.59 ns, 
where the color image indicates the temperature distribution and the arrows indicate the heat flux. 
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wide wavelength shift of 42% for double FRs is observed in the M-IR 
region. This tunable effect is due to the phase transition from 
amorphous to crystalline. A heat transfer model is built to resolve 
the transient temperature variation in the structure during a photo - 
thermal process. Our model predicts that a fast phase transition time 
of 0.36 ns can be reached under a weak pump intensity of 9.6 /nW/ 
^m 2 .This work introduces a new approach for the study of ultralow- 
power and ultrafast photonic nanodevices. It may open up avenues 
for designing active multiple -wavelength slow light devices, plasmo- 
nics switching and sensors. 

Methods 

For the numerical calculations, we used the Finite Element Method by means of 
COMSOL Multiphysics. The unit cell of the phase-change metamaterial is shown in 
Fig. 1 . The incident light is normal to the x-y plane with E field polarized in y direction. 
To account for the periodic nature of the metamaterial, the model boundary at 

x= + -y and y = ±L y is set to condition of perfect magnetic conductor and perfect 

electric conductor, respectively for normal incidence. Scattering boundary conditions 
are applied in the z direction. The optical constant of gold is taken from the data of 
Johnson and Christy 46 .The dielectric function of Ge 2 Sb 2 Te 5 in different phases are 
obtained from the experimental data in 42 . 
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